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Abstract 

We derive a simple analytical scaling function which embodies the predictions of high 
density QCD on the energy, centrality, rapidity, and atomic number dependences of hadron 
multiplicities in nuclear collisions. Both centrality and rapidity dependences of hadron 
Q ' multiplicity in Au — Au collisions as measured at RHIC at y/s = 130 GeV are well 

described in this approach. The centrality and rapidity dependences of hadron multiplicity 
at \fs = 200 GeV run at RHIC are predicted; the variation of these dependences with 
.^J \ energy appear different from other approaches, and can be used as an important test of 

^ ■ the ideas based on parton saturation and classical Chromo-Dynamics. 

Relativistic heavy ion collisions allow to probe QCD in the non-linear regime of high parton 
density and high color field strength. It has been conjectured long time ago that the dynamics 
of QCD in the high density domain may become qualitatively different: in parton language, 
this is best described in terms of parton saturation fl], |2|, ||| , and in the language of color fields 
- in terms of the classical Chromo-Dynamics [|], |5], |], 0. In this high density regime, the 
transition amplitudes are dominated not by quantum fluctuations, but by the configurations of 
classical field containing large, ~ l/ct s , numbers of gluons. One thus uncovers new non-linear 
features of QCD, which cannot be investigated in the more traditional applications based on 
the perturbative approach. The classical color fields in the initial nuclei (the "color glass con- 
densate" 0, can be thought of as either perturbatively generated, or as being a topologically 
non-trivial superposition of the Weizsacker- Williams radiation and the quasi-classical vacuum 
fields 

Recently, with the advent of RHIC, the heavy ion program has entered the era of collider 



experiments. First RHIC results on hadron multiplicities have been presented []10 , pd] , |12j , [1 
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The established experimental centrality dependence is so far in accord with predictions pj based 
on high density QCD. (For alternative approaches, see |14|, [15|, [L7|]). Rapidity dependence 



of the hadron multiplicity has recently been presented as well, and it is important to check 
if it agrees with the high density QCD calculations as well. It is also most timely to make 
predictions for the higher energy, a/s = 200 GeV run. These are the objectives of this paper. 

Let us begin by presenting a brief and elementary introduction into the concept of parton 
saturation [|l|. Consider an external probe J interacting with the nuclear target of atomic 
number A (see Fig. 1). At small values of Bjorken x, by uncertainty principle the interaction 
develops over large longitudinal distances z ~ l/(mx), where m is the nucleon mass. As soon as 
z becomes larger than the nuclear diameter, the probe cannot distinguish between the nucleons 
located on the front and back edges of the nucleus, and all partons within the transverse area 
~ 1/Q 2 determined by the momentum transfer Q participate in the interaction coherently. The 
density of partons in the transverse plane is given by 

^ xG A (xQ 2 ) ^ Al/S 
ttR a 

where we have assumed that the nuclear gluon distribution scales with the number of nucleons 
A. The probe interacts with partons with cross section a ~ a s /Q 2 'i therefore, depending on the 
magnitude of momentum transfer Q, atomic number A, and the value of Bjorken x, one may 
encounter two regimes: 

• ap A C 1 - this is a familiar "dilute" regime of incoherent interactions, which is well 
described by the methods of perturbative QCD; 

• op A 3> 1 - in this regime, we deal with a dense parton system. Not only do the "leading 
twist" expressions become inadequate, but also the expansion in higher twists, i.e. in 
multi-parton correlations, breaks down here. 

The border between the two regimes can be found from the condition ap A — 1; it determines 
the critical value of the momentum transfer ( "saturation scale" [|T], |J) at which the parton system 
becomes to look dense to the probeQ: 

2 xG A (x,Q 2 ) 

Qs ~ a s —52 • ( 2 ) 

ttK a 

In this regime, the number of gluons from (fj) is given by 

xG A (x,Cfi) ~ -L- C£& A , (3) 

where Q 2 S R A ~ A. One can see that the number of gluons is proportional to the inverse of 
a s (Q 2 ), and becomes large in the weak coupling regime. In this regime, as we shall now discuss, 
the dynamics is likely to become essentially classical. 



1 Note that since Q 2 S ~ A 1 / 3 , this expression in the target rest frame can also be understood as describing a 
broadening of the transverse momentum resulting from the multiple re-scattering of the probe. 
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Figure 1: Hard probe inter- 
acting with the nuclear target 
resolves the transverse area ~ 
1/Q 2 (Q 2 is the square of the 
momentum transfer) and, in 
the target rest frame, the lon- 
gitudinal distance ~ l/(mx) 
[m is the nucleon mass and x 
- Bjorken variable). 



The condition (0) can be re-derived in a way |18| which illustrates the link between the 
saturation and classical Yang-Mills fields ||. As a first step, let us re-scale the gluon fields in 
the Lagrangian 

C = -\d%G% + £ tfiil.D, - m f )q a f] (4) 
4 / 



as follows: A a u 



Al = gA°. In terms of new fields, G%, = gG^ = d^A a u - d v A* + f^A^Al, 
and the dependence of the action corresponding to the Lagrangian (|j) on the coupling constant 
is given by 

S ~ I j 2 G%G% d 4 x. (5) 



Let us now consider a classical configuration of gluon fields; by definition, G a in such a con- 
figuration does not depend on the coupling, and the action is large, S 3> h. The number of 
quanta in such a configuration is then 



a.. 



(6) 



where we re-wrote (|5]) as a product of four-dimensional action density p 4 and the four 
dimensional volume V4. 

The effects of non-linear interactions among the gluons become important when d^A^ ~ 
(this condition can be made explicitly gauge invariant if we derive it from the expansion of a 
correlation function of gauge-invariant gluon operators, e.g., G 2 ). In momentum space, this 
equality corresponds to 

A* ~ (G 2 ) 1/2 = VpI; (7) 



A 2 



Q s is the typical value of the gluon momentum below which the interactions become essentially 
non-linear. 

Consider now a nucleus A boosted to a high momentum. By uncertainty principle, the gluons 
with transverse momentum Q s are extended in the longitudinal and proper time directions by 
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~ 1/Qs', since the transverse area is ttR\, the four-volume is V4 ~ TcRj^/Q 2 . The resulting 
four-density from @ is then 

p 4 ~ a s -ft ~ a s —-2- ~ Q s , (8) 

where at the last stage we have used the non-linearity condition (0), p 4 ~ Qj. It is easy to 
see that coincides with the saturation condition (0), since the number of gluons in the 
infinite momentum frame N g ~ xG(x, Q 2 )- This simple derivation illustrates that the physics 
in the high-density regime can potentially be understood in terms of classical gluon fields. This 
correspondence allowed to formulate an effective quasi-classical theory j4|, which is a subject 
of vigorous investigations at present (see, e.g., 0, 0). 

In nuclear collisions, the saturation scale becomes a function of centrality; a generic feature 
of the quasi-classical approach - the proportionality of the number of gluons to the inverse of 
the coupling constant @ - thus leads to definite predictions on the centrality dependence 
of multiplicity: 

dN 

— = c N part xG(x,Q 2 ), (9) 

where xG(x, Q 2 ) ~ \/a s {Q 2 s ) ~ \n(Q 2 / Aq CD ), N part is the number of participants, and c is the 
"parton liberation" coefficient, which was determined from the experimental multiplicity to 
be 

c= 1.23 ±0.20. (10) 



This value is on the order of unity fll9"| , and agrees with the lattice [50 and analytical pi] 
calculations. The prediction (§) so far is in accord with the data coming from RHIC [|K], [LlJ 



12], 13"]. Let us now make predictions for y/s = 200 GeV collisions based on this picture. 



The energy dependence of the hadron production is determined by the variation of saturation 
scale Q s with Bjorken x = Q s / yfs. This variation, in turn, is determined by the x— dependence 
of the gluon structure function. In the saturation approach, the gluon distribution is related to 
the saturation scale by Eq.(0). A good description of HERA data is obtained with saturation 
scale Q 2 S = 1 4- 2 GeV 2 with W - dependence (W = y/s is the center-of-mass energy available 
in the photon-nucleon system) 

Ql oc W\ (11) 

where A ~ 0.25 -4- 0.3. In spite of significant uncertainties in the determination of the gluon 
structure functions, perhaps even more important is the observation JEJ that the HERA data 
exhibit scaling when plotted as a function of variable 



91 

Ql \x Q ) 



(12) 



where the value of A is again within the limits A ~ 0.254-0.3. In high density QCD, this scaling 
is a consequence of the existence of dimensionful scale |l], J3J) 

Q 2 s (x)=Ql(x /x) x . (13) 
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The formula (|) can be equivalently re-written as (see J7| for details) 



dN q Q * iu) 
*7 Sa ^M)' (14) 

where Sa is the nuclear overlap area, determined by atomic number and the centrality of 
collision. 

The energy dependence of multiplicity is thus (up to a logarithmic correction, which is small 
for the energy interpolation that we make) given by 

^(V5 = 200GeV)~(^T) ^(v^= 130 GeV) = (1.10-1.14) ^(v^= 130 GeV), (15) 
dr] Viol)/ ar] dr] 

indicating that the hadron multiplicity is expected to raise in the y/s = 200 GeV run by about 
10 -4- 14%5 Taking the PHOBOS number for charged multiplicity at y/s = 130 GeV for the 6% 
centrality cut, dN/dr] = 555 ± 12(stat) ± 35(syst) f25j, one gets a prediction of 

dN/dr] = 616 4- 634 (16) 

for the y/s = 200 GeV (these are the central values not taking into account the error bars of 
the y/s = 130 GeV measurement). 

One can also try to extract the value of A from the energy dependence of hadron multiplicity 
measured by PHOBOS at y/s = 130 GeV and at at y/s = 56 GeV; this procedure yields 
A ~ 0.37, which is larger than the value inferred from the HERA data (and is very close to the 
value A ~ 0.38, resulting from the final-state saturation calculations |L5|1 ). 

Let us now proceed to the calculation of the (pseudo)rapidity and centrality dependences. 
We need to evaluate the leading tree diagram describing emission of gluons on the classical level, 
see Fig. 2. Let us introduce the unintegrated gluon distribution <pa(%, k% ) which describes the 
probability to find a gluon with a given x and transverse momentum kt inside the nucleus A. 
As follows from this definition, the unintegrated distribution is related to the gluon structure 
function by 

2 

fPt 

xG A (x,p 2 t )=J dk\ <p A (x,kt)] (17) 

when pj > Q 2 S , the unintegrated distribution corresponding to the bremsstrahlung radiation 
spectrum is 

a., 1 ,„„x 



<p A (x,k*) ~ ^ ^ 



In the saturation region, the gluon structure function is given by (|3[); the corresponding unin- 
tegrated gluon distribution has only logarithmic dependence on the transverse momentum: 

Pa(x, kf) ~ — ; k\<Ql (19) 



2 At the time when this paper is being finalized, the result of the PHOBOS multiplicity measurement for 
central collisions has just been announced - see Note added. 
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Figure 2: The Mueller dia- 
gram for the classical gluon ra- 
diation. 



V x 2'(P t - k t) > 



where Sa is the nuclear overlap area, determined by the atomic numbers of the colliding nuclei 
and by centrality of the collision. 

The differential cross section of gluon production in a AA collision can now be written down 
as § H 

do AttN c 1 f 2 , 3 . . , 2 

d^) = N 2 -l p 2 t J 1 a& VA\X\,h) (p A (x2,(p-k) t ), (20) 

where x± } 2 = (pt/y/s) exp(±^), with rj the (pseudo)rapidity of the produced gluon; the running 
coupling ot s has to be evaluated at the scale Q 2 = max{k 2 , (p — k) 2 }. The rapidity density is 
then evaluated from (p0|) according to 



dN 



dy oaa 




(21) 



where oaa is the inelastic cross section of nucleus-nucleus interaction. 

Since the rapidity y and Bjorken variable are related by In 1/x = y, the x— dependence of 
the gluon structure function translates into the following dependence of the saturation scale Q 2 
on rapidity: 

Q%s; ±y) = Q 2 s {s- y = 0) exp(±Ay). (22) 

As it follows from (^) , the increase of rapidity at a fixed W = \fs moves the wave function 
of one of the colliding nuclei deeper into the saturation region, while leading to a smaller gluon 
density in the other, which as a result can be pushed out of the saturation domain. Therefore, 
depending on the value of rapidity, the integration over the transverse momentum in Eqs. 
(^I]),(^T|) can be split in two regions: i) the region Aqcd < h < Q s ,min in which the wave 
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functions are both in the saturation domain; and ii) the region A << Q s , m in < h < Q s ,max 
in which the wave function of one of the nuclei is in the saturation region and the other 
one is not. Of course, there is also the region of k t > Q s ,max, which is governed by the usual 
perturbative dynamics, but our assumption here is that the role of these genuine hard processes 
in the bulk of gluon production is relatively small; in the saturation scenario, these processes 
represent quantum fluctuations above the classical background. It is worth commenting that 
in the conventional mini-jet picture, this classical background is absent, and the multi-particle 
production is dominated by perturbative processes. This is the main physical difference between 
the two approaches; for the production of particles with p t >> Q s they lead to identical results. 

To perform the calculation according to (pl"|),(f2"0|) away from y = we need also to specify 
the behavior of the gluon structure function at large Bjorken x (and out of the saturation 
region). At x — > 1, this behavior is governed by the QCD counting rules xG(x) ~ (1 — x) 4 , 
so we adopt the following conventional form: xG(x) ~ x~ x (1 — x) 4 . 

We now have everything at hand to perform the integration over transverse momentum in 
(pl|), (p0|); the result is the following: 



= const Sa Q 

dy 



2 

s,min 



In 



'Q 



2 

s,min 



A 2 

^QCD 



1 + - In 
2 



'Q 



i 

s,max 



Q 



2 

s,min , 



1 



Q 



S.IIHIX 



(23) 



where the constant is energy-independent, Sa is the nuclear overlap area, Q 2 = Q 2 s (s; y = 0), 



and Q 



s,min(max) 



are defined as the smaller (larger) values of (P2|); at y — 0, Q 



) 2 . 

' s,min 



Q 



2 

s,max 



Ql( s ) = Qs( s o) ( s / s o) X ^ 2 - The first term in the brackets in (p3|) originates from the region in 
which both nuclear wave functions are in the saturation regime; this corresponds to the familiar 
~ (l/cts) Q 2 S R\ term in the gluon multiplicity, cf. Eq.(|3]). The second term comes from the 
region in which only one of the wave functions is in the saturation region. The coefficient 1/2 
in front of the second term in square brackets comes from kt ordering of gluon momenta in 
evaluation of the integral of Eq . (^0|) . 

The formula (^) has been derived using the form ([H]) for the unintegrated gluon distri- 
butions. We have checked numerically that the use of more sophisticated functional form of 
if a taken from the saturation model of Golec-Biernat and Wusthoff |22|] in Eq.(|20D affects the 
results only at the level of about 3%. 



Since SaQ 2 ~ N part (recall that Q 2 S ^> Aq CD is defined as the density of partons in the 



transverse plane, which is proportional to the density of participants), we can re-write 
the following final form 



3 in 



dN 
dy 



C N part I — 



-AM 



hi 



Qt 



A 2 

s QCD , 



(24) 



with Q 2 (s) = Q 2 (s ) (s/so) x / 2 . This formula is the central result of our paper; it expresses 
the predictions of high density QCD for the energy, centrality, rapidity, and atomic number 
dependences of hadron multiplicities in nuclear collisions in terms of a single scaling function. 
Once the energy-independent constant c ~ 1 and Q 2 (so) are determined at some energy so, Eq. 
(ID) contains no free parameters. (The value of A, describing the growth of the gluon structure 
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functions at small x can be determined in deep-inelastic scattering; the HERA data are fitted 
with A ~ 0.25 -i- 0.3 [p2fl ). At y = the expression (|2~3"D coincides exactly with the one derived 
in j7|, and extends it to describe the rapidity and energy dependences. 

Before we can compare (^3|) to the data, we have to take account of the difference between 
rapidity y and the measured pseudo-rapidity r\. This is done by multiplying ( p3|) by the Jacobian 
of the y 7] transformation; a simple calculation yields 



h{rj;p t ;m) 



cosily 



^f-*- + smh rj 



(25) 



pi 



where m is the typical mass of the produced particle, and p t is its typical transverse momentum. 
Of course, to plot the distribution (|2~4]) as a function of pseudo-rapidity, one also has to express 
rapidity y in terms of pseudo-rapidity rj; this relation is given by 



y(ij;pt;m) = - In 



m \ Vt + sinh 2 rj + sinh rj 

Pt 



Pt 



+ sinh r] — sinh r\ 



(26) 



obviously, h(r];p t ;m) = dy(r)',pt',m)/dr]. 

We now have to make an assumption about the typical invariant mass m of the gluon 
mini-jet. Let us estimate it by assuming that the slowest hadron in the mini-jet decay is the 
p-resonance, with energy E p = [m 2 + p 2 t + p 2 p J 1 / 2 , where the z axis is pointing along the 
mini-jet momentum. Let us also denote by Xi the fractions of the gluon energy q carried by 
other, fast, i particles in the mini-jet decay. Since the sum of transverse (with respect to the 
mini-jet axis) momenta of mini-jet decay products is equal to zero, the mini-jet invariant mass 
m is given by 



rn 



jet 



= m 



; Xiq + Ep) 2 - X i<lz + Ppr 



2 V Xiq z 



[171 



p.t 



Pp,z) — ^Qs 



m. 



effi 



(27) 



where m 



{ml + p 2 



vl/2 



1 and q ~ q z = Q s . Taking 



. In Eq. (pTD we used that J2i x i 
P P ,z ~ P P ,t ~ 300 MeV and p mass, we obtain m e // ~ 0.5 GeV. 

We thus use the mass m 2 ~ 2Q s m e ff ~ Q s ■ 1 GeV in Eqs. ( p5] , p6] ) . Since the typical 
transverse momentum of the produced gluon mini-jet is Q s , we take p t = Q s in (p5|). The 
effect of the transformation from rapidity to pseudo-rapidity is the decrease of multiplicity at 
small rj by about 25 — 30%, leading to the appearance of the ~ 10% dip in the pseudo-rapidity 
distribution in the vicinity of rj = 0. We have checked that the change in the value of the 
mini-jet mass by two times affects the Jacobian at central pseudo-rapidity to about ~ 10%, 
leading to ~ 3% effect on the final result. 

The results for the Au — Au collisions at ^J~s =130 GeV are presented in Figs § and |j. In 
the calculation, we use the results on the dependence of saturation scale on the mean number 
of participants at y/s = 130 GeV from 0, see Table 2 of that paper. The mean number of 
participants in a given centrality cut is taken from the PHOBOS paper ||25||. One can see that 
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Figure 3: Centrality dependence of charged hadron production per participant at different 



pseudorapidity r\ intervals in Au — Au collisions at ^fs = 130 GeV; the data are from [BS . 
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Figure 4: Pseudo-rapidity dependence of charged hadron production at different cuts on 
centrality in Au — Au collisions at y^s = 130 GeV; the data are from [^5 . 
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Figure 5: Centrality dependence of charged hadron production per participant at different 
pseudorapidity rj intervals in Au — Au collisions at y/s = 200 GeV. 
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Figure 6: Pseudo-rapidity dependence of charged hadron production at different cuts on cen- 
trality in Au — Au collisions at = 200 GeV. 
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both the centrality dependence and the rapidity dependence of the \/s = 130 GeV PHOBOS 
data are well reproduced below r\ ~ ±4. The rapidity dependence has been evaluated with 
A = 0.25, which is within the range A = 0.25 -j- 0.3 inferred from the HERA data p2"| . The 
discrepancy above r\ ~ ±4 is not surprising since our approach does not properly take into 
account multi-parton correlations which are important in the fragmentation region. 

Our predictions for Au — Au collisions at = 200 GeV are presented in Figs. |5] and || 
The only parameter which governs the energy dependence is the exponent A, which we assume 
to be A ~ 0.25 as inferred from the HERA data. The absolute prediction for the multiplicity, as 
explained above, bears some uncertainty, but there is a definite feature of our scenario which is 
distinct from other approaches. It is the dependence of multiplicity on centrality, which around 
7] = is determined solely by the running of the QCD strong coupling 0. As a result, the 
centrality dependence at y/s = 200 GeV is somewhat less steep than at y/s = 130. While the 
difference in the shape at these two energies is quite small, in the perturbative mini-jet picture 
this slope should increase, reflecting the growth of the mini-jet cross section with energy. 
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saturation (solid line) and 
"soft plus hard" (dashed line) 
calculations of the centrality 
dependence of charged hadron 
multiplicity at \fs = 200 GeV; 
also shown is the "soft plus 
hard" calculation at y/s = 130 
GeV (lower dashed curve). . 
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Let us confirm this statement quantitatively. The following simple parameterization J7| was 



found to describe the data [10], O, [12[ quite well: 



dN 
drj 



n 



VP 



(Neon), 



(2f 



where (N part ( coU )) are the average numbers of participants (collisions) for a given centrality, 
and X(s) is an s— dependent parameter reflecting the relative strength of "soft" and "hard" 
components in multi-particle production. At yfs = 130 GeV, the value deducted from the 
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central multiplicity was X(130 GeV) = 0.09 ±0.03, and at = 56 GeV, it was X(56 GeV) = 
0.05 ± 0.03. In the perturbative picture, the "hard" component of multiplicity is proportional 
to the mini-jet production cross section. At rj — 0, the energy dependence of the mini-jet 
cross section is determined by the square of the gluon structure function xG(x) (note that 
in the saturation approach, as we discussed above, it is the first power of xG(x)). Taking 
xG(x) ~ x~ A , with A ~ 0.3, translates then in the increase of the mini-jet production cross 
section by factor of (130/56) 2A ~ 1.7 in going from y/s = 56 GeV to y/s = 130 GeV, in accord 
with the ratio of the extracted values of X(s) quoted above. Extrapolating the same energy 
dependence to y/s = 200 GeV, we estimate 

X(y/s = 200 GeV) = ( X(y/s = 130 GeV) ~ 0.12 ± 0.02, (29) 

V 130 / 

where the error bar reflects only the uncertainty in the energy extrapolation of the central 
values of X. The value (|29|), together with Eq.([28"D, represents the prediction of the "soft plus 



hard" approach for the centrality dependence at y/s = 200 GeV. Using the value (^) together 



with the average numbers of participants from |7| (see |26| for details of the Glauber model 
we use) (N part ) = 339 and collisions (N co u) = 1049 in the 6% centrality cut at y/s = 200 
GeV, and the pp multiplicity of n pp = 2.43 which follows from the parameterization of the data 
n pp = 2.5 - 0.25 ln(s) + 0.023 ln 2 (s) J27|, [2jJ, we get the central value 

dN , 

—On = ) ~ 668, 30 
dr\ 

which represents a ~ 20% increase compared to y/s = 130 GeV. This value is larger than the 
one predicted above on the basis of the initial-state saturation approach; however the main 
difference between the two scenarios will be in the centrality and rapidity dependences of hadron 
multiplicity. 

To estimate the difference in the shape of centrality dependence, let us consider the case 
when the multiplicities in central collisions are forced to be the same and equal to dN/df](r] = 
0) ~ 634, which is the upper limit of our prediction (|16"D and below the central value (|30|) (note 
that this represents the "worst case" in which the value of multiplicity in central collisions 
alone cannot be used to distinguish the two approaches). In the "soft plus hard" model, this 
multiplicity would then correspond to X = 0.106. In Fig. we compare the resulting centrality 
dependences in the saturation and in the "soft plus hard" calculations. While the difference 
between the two curves is not large, the variation of the predictions with energy and rapidity 
could help to discriminate between the two approaches; for comparison, we also show the result 
of the "soft plus hard" calculation at y/s =130 GeV. 

To summarize, we have extended the previous analysis of multi-particle production in the 
quasi-classical picture to include the pseudo-rapidity dependence, and found good agreement 
with RHIC data at y/s = 130 GeV. If this agreement persists at higher energies, one may 
conclude that the dynamics of quasi-classical color fields, encoded in our approach in terms 
of a simple scaling function (0), describes well the gross features of multi-particle production 
already at RHIC energy. The experimental study of centrality and rapidity dependences at 
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y/s = 200 GeV will thus be a very important step in establishing the presence of high density 
QCD effects in relativistic nuclear collisions. The clear discrimination between this and other 
approaches may however require the study of more "microscopic" observables. 

Note added: While we were finalizing this paper, the first y/s = 200 GeV multiplicity mea- 
surement by the PHOBOS Collaboration had been announced p8 . The ratio of multiplicities 



in the central (0 — 6% centrality cut) events at y/s = 200 and y/s = 130 GeV has been reported: 
-R(200/130) = 1.14 ± 0.06. While this is in accord with the prediction based on the saturation 
approach given above, _R(200/130) = 1.10 ~ 1.14, it is also consistent with the 20% increase 
which we expect in the "soft plus hard" scenario. (We did not attempt to adjust to the central 
value of the announced data our calculations shown in Figs |5] and |] and reported (e.g., f29fl ) 
before the multiplicity in head-on collisions had been available). The crucial test of the quasi- 
classical approach will thus be provided by the centrality and pseudo-rapidity dependences of 
hadron multiplicities. 
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98CH10886). The research of E.L. was supported in part by the Israel Science Foundation, 
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